The ground state energy of a low density Bose gas: a second order upper bound 
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Consider N bosons in a finite box A = [0,L] 3 C R 3 interacting via a two-body nonnegative soft 
potential V = XV with V fixed and A > small. We will take the limit L, N — * oo by keeping the 
density g = N/ L 3 fixed and small. We construct a variational state which gives an upper bound on 
the ground state energy per particle e 
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1 < S x < 1 + C\ 



Here a is the scattering length of V and thus depends on A. In comparison, the prediction by 
Lee- Yang and Lee-Huang- Yang fjl] asserts that S\ = 1 independent of A. 
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INTRODUCTION 

Although Bose-Einstein condensation has been firmly 
established since the experiments 0, rigorous un- 
derstanding of the Bose gas starting from the many- 
body Schrodingcr equation is still in a very rudimentary 
stage and many theoretical predictions at present are still 
based on uncontrolled approximations. Notable excep- 
tions are the rigorous derivations of the Gross-Pitaevskii 
equation in both the stationary and the dynamical set- 
tings [!, HH- In particular, in the limit of low density, 
the leading term of the ground state energy per particle 
of an interacting Bose gas was proved by Dyson (upper 
bound) d and Lieb-Yngvason (lower bound) [13| to be 
Airag where a is the scattering length of the two-body 
potential and g is the density. The famous second order 
correction to this leading term was first computed by 
Lee- Yang [l(| (see also Lee-Huang- Yang Q and the re- 
cent paper by Yang 16] for results in other dimensions). 
In this paper, we construct a trial function with a second 
order term in the energy which, up to a constant factor, 
is the same as predicted in To present this result, 

we now introduce our setup rigorously. 

Consider N interacting bosons in a finite box A = 
[0, L] 3 C R 3 with periodic boundary conditions. Let V 
be a smooth, radially symmetric, nonnegative potential 
with fast decay. The two-body interaction V is given by 
V = XV with A a small constant. We will first take the 
limit L, N — > oo by keeping the density g = N/L 3 fixed. 
In the limit g — > 0, the leading term for the ground state 
energy per particle of this system is Airga, where a is the 
scattering length of the potential V. The Lee- Yang's 
prediction of the energy per particle up to the second 
order term is given by 

128 =(e« 8 ) 1/a + -- 
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potential approximation 0, Q and the "binary collision 
expansion method" 0. One can also obtain fl]) by per- 
forming the Bogoliubov Q approximation and then re- 
placing the integral of the potential by its scattering 
length [8|1 , Another derivation of |T]) was later given by 
Lieb till using a self-consistent closure assumption for 
the hierarchy of correlation functions. Although these 
approaches gave the same answer for the second order 
term ((T|), it is nevertheless difficult to extract rigorous 
results on the energy using these ideas. In fact, the first 
rigorous upper bound on the energy to the leading order 
by Dyson [J] was based on a completely different con- 
struction. The proof of Lieb-Yngvason [13| on the lower 
bound of the energy to the leading order was also very 
different from the earlier approaches. 

The trial wave function of Dyson Q also shows that 
the next order correction in energy for hard core bosons 
is bounded from above by C(ga 3 ) 1 ^ 3 . The same upper 
bound for soft potentials was obtained in [T^]. In this 
paper, we construct a variational state which gives a 
rigorous upper bound on the ground state energy per 
particle 



e < Anga 
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1+ (ga 3 )^S x +0(g 2 \\ogg\) (2) 



The approach by Lee- Yang [lOj is based on the pseudo- 



with S\ < 1 + CX. The second order term in this upper 
bound ([2|) is of the same form as the conjectured one 
|T|), up to a positive correction in the constant of order 
A. This constant C and the constant in the error term in 
depends on the details of the interaction potential, 
in particular our proof uses that V is a soft potential. 

The trial state in this paper consists only of conden- 
sate particles and of non-condensate particle pairs with 
momenta k,—k, reminiscent of the original idea of Bo- 
goliubov. In our computation, however, the interactions 
between non-condensate particle pairs are also relevant. 
Our trial state does not have a fixed number of particles, 
but it is a state in the Fock space with expected number 



2 



of particles equal to N. It is similar to the trial states 
used by Solovej [1 51 ] to give rigorous upper bounds to 
the ground state energies of the one and two-component 
charged Bose gases in the high density limit. 

Variational trial states with particle pairs have been 
used earlier in the context of the low density Bose gas by 
Girardeau and Arnowitt 0]. Their state, however, had a 
fixed number of particles which slightly complicated the 
calculation (the details are available only in the unpub- 
lished Ph.D. dissertation of Girardeau). The variational 
formula we obtain is nevertheless the same as theirs up 
to lower order terms due to the choice of a different en- 
semble. However, in [6j, the solution of the minimiza- 
tion problem was given only implicitly as a solution to 
a nonlinear integral equation and thus the energy was 
not evaluated explicitly. In our work, we identify the 
presumed main terms from the calculations of each indi- 
vidual terms in the energy. This enables us to find the 
minimizer for the main terms of the energy. By choosing 
the minimizer of the main part as our trial state, we a- 
posteriori justify that the neglected terms are indeed of 
lower order and thus giving a rigorous upper bound on 
the energy. We believe that the difference between the 
energy of our state and that of the true minimizer of the 
full functional is of lower order. 



SETUP 

We work in a finite box A = A L = [0, L} 3 C R 3 with 
periodic boundary conditions. Its dual space is A* := 
(^Z) 3 . For a continuous function F on IR 3 , we have 

lim -L Y F{p) = Hm-^V F(p) 



peA* 



d 3 p 
(27)3 



pSA* 



F(p) 



The Fourier transform of an arbitrary function f(x) on 
A is defined as 

1 1 peA* 



fp = I er^f(x)dx, 

I A 



Note that the Fourier transform depends on A, a fact 
that is omitted from the notation. In most cases we 
will take Fourier transforms of sufficiently decaying func- 
tions, so that their A dependence is negligible in the limit 
L — > oo. Since V(x) is real and symmetric, we have that 
V u is real and 

V, = • 



indicates whether it is the momentum or position space 
delta function. We also simplify the notation 

E-E 

p peA' 

i.e. momentum summation is always on the whole A*. 

Notice that the choice of the Sa* ensures that in the 
L — > oo limit, this delta function converges to the usual 
continuum delta function 8{p) in momentum space with 
respect to the measure d 3 p/(27r) 3 : 

fern Y.np)^{p- q ) 

1 1 pGA* 



d 3 p 

(27T) 5 



F(p)S(p -q) = F(q) 



(where 5(p) is defined by the last formula). 

Using the formalism of second quantization, we work 
on the bosonic Fock space T — ©^Lo^® 3 ™ built upon 
the single particle Hilbert space H = l 2 {A*), where Tl®° n 
denotes the symmetric tensor product of n copies of Ti.. 
The vacuum is denoted by |0). We consider bosonic an- 
nihilation and creation operators, a,k, , for any k G A*, 
with the usual canonical commutation relations (CCR): 

[op, 5+] = ~a p a+ - a+a p = 8(p-q)=^ ftLTwle. 
The Hamiltonian of the system is given by 
H = W\ E p2fi p h P + 5lX|3 E V u a+a+a p - u a q+u . 



2\kf 



where the first term is the kinetic energy, the second 
term is the interaction energy of the particles in appro- 
priate physical units. It is more convenient to redefine 
the bosonic operators as 



1 



1 



k ' 



i.e., from now on we assume that 

r +i + + ( 1 i£p = q 

1 p i 1 pi ip o otherwise. 



Thus the Hamiltonian is given by 
H = Y,P 2a t a P + 2jX| E ^ a P 



" q ""p—u^q+u 



",: ",. „«., i ,. ■ (3) 



We also have 



|Af E ' 

1 1 peA* 



,tp-x 



8m? (x), 



e l P x dx = 5A*{p) 



where 5r3 (x) is the usual continuum delta function and 
<5a* is the lattice delta function, i.e. 5a* (p) = |A| = L 3 
if p = 0, and 8 A * (p) = if p £ A* \ {0}. We will ne- 
glect the subscripts in the delta functions, the argument 



THE TRIAL STATE 

Let Ck denote a family of complex numbers 
parametrized by k G A* \ {0} with the property that 
led < 1 and Ck = c-k- We define a state 



> 2 E* 



(4) 
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where iVo is a positive real number. The parameters Ck 
and No will be fixed later on. 

Fix a small positive number g which will be the density 
of the system and let N denote 

N:=q\A\ . 

Define the expectation w.r.t. "J by 



where (•,•) denote the standard L 2 inner product. We 
shall fix the parameters in \I> so that the expected num- 
ber of particles is given by N 



N = ( a ™ a r* 



(5) 



m£A* 



Let E = (i?)* be the energy. 

Let 1-wbe the zero energy scattering solution to the 
potential V 



A(l-w) + ~V(l-w) =0 



(6) 



on R 3 with < w < 1 and w(x) — ► as \x\ — > oo. Then 
the scattering length is defined by 



a:=^- I \'(.r)(l - ,rU)Hl.r 



—— lim / V(x)(l — w(x))dx . 
on l^oc J A 



(7) 



ft is well-known that 

8ira < { V(x)dx = lim / V(x)dx = lim V , (8) 

J R3 L->caJ A L^oo 

where, in the last step, we recall that the definition of 
the Fourier transform depends on L. 
Define the number h by 



Theorem 1 Let V(x) > 0, V ^ 0, be a non-negative 
radially symmetric smooth function with a decay V(x) < 
C(l + \x\)- 3 - 5 for some S > 0, and set V(x) = XV (x). 
Then for X small enough, we have, in the limit g — > 0, 
the following estimate 



E = AngNa + Q + 0(Ng 2 \ log g\) 



(12) 



for the energy of the trial state @ with an appropriate 
choice of Ck and No, under the constraint ([5]). Here 
Q = Q(h) is given by 



-$(M(a 3 g) 1/J 

7T 



Q(h) = AiraNg 



and the constant in the error term in (|12p depends on X 
and V. 

The assumptions on V can certainly be relaxed but we 
do not aim at identifying the optimal conditions. 
A direct calculation gives 



*(0) 

thus at h = we obtain 

Q(0) = inaNg- 



/5L2 
T5~ 
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(a 3 g) 



1/2 



Moreover, a simple calculation also shows that 

2y 1/2 dy 



f 

< $'(0) = / 
Jo 



< oo, 



/o (y + lWv(y + 2) 

thus infinitesimally $(0) < if < h -€i 1 and 

Q{h) = Q(0) + 0(h) = 0(0) + O(A) 
by pTJj) . In fact, it is also easy to see that 
$(0) < 

holds for any /i > 0. Thus our trial state delivers a 
second order term with an explicit constant that is bigger 
than the Lee- Yang prediction 0, [l(| by a factor (1 + 
O(A)) for small coupling constant A. 



h = 3__ 1 

8ira 



(9) 



from © it follows that h > if L is sufficiently large. 
Recall that V = XV with V" being fixed. The scattering 
length a can be computed via the Born series for small 
A. In particular, we will show in Lemma [3] that h is of 
order A 



h = O(A) 



(10) 



Define the function 

*(h) := / &yy 1/2 
Jo 



x (V (y + 2h)(y + 2 + 2h) - (y + 1 + 2h) + ^) . 

ai) 

One can check that this integral is convergent for h > 0. 
Our main result is the following theorem. 



COMPUTATION OF THE ENERGY 

We start the proof of the main theorem by the follow- 
ing Lemma. We first define a few quantities: 




(13) 



Op | |L. r | 



-|A| 



(l-| Cp |2)(l-| Cr P) 

1 I/ |c P | 2 (l + 3|c p | 2 ) + V 2p \c p \ 2 (l + \c p \ 2 ) 



2IAI 



(l-|c P | 2 ) 2 



(14) 
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Lemma 2 The energy E — (if)* of the state @ under 
the constraint ([5]) is given by E = Em + CI2 + ^4, where 



(15) 



+ 



E 

1 



(p 2 + gV P ) 



i-M 2 



- i?ec„ 



2|A| ? ^ _r fl 



r#0,±p 



£. \ " 

|A| ^ 1 

1 1 m^O 



I I 



(l-| Cp P)(l-|c r P) 



Rec r , 



1 - c 



Remark: We shall see later on that the first four terms 
in the main term, Em, are of order Ng, the fifth one is 
of order Ng 3 / 2 . Each term in the error terms Q2 and ^4 
is at most of order Ng 2 \ log g\ except the last term in ft 4 
which is non-extensive. 

Proof. We collect a few elementary facts from direct 
calculations. Similar formulas in an abstract setup ap- 
plicable in general, not only for the translation invariant 
case, were presented in (l5j . 



,c k a k a_ k c k a k a_ k 



10) = r 



1 



_ o N 



\ c k\ 

n 



for fc^O, 
1 



k^o (l ^ |cfe| : 



,1/2 ' 



(a£a ) 9 = N ■ 
(a+a+a a )<p = iV 2 , 

and for m/0, 

(ftmOm)* = 



(a m a_ m )y — (a m a_ m )<s> 



l-\c m ? 

(-111 



* l-|c m | 2 

( a ro a m a -ro a -m)* = 

|c m | 2 (l + |c m | 2 ) 



(a m a_ m a m a_ m )* — 



l-|c m | 2 ) 2 
21c 14 



Moreover, 



\ a m a m a m a m) * — (1 — |c | 2 ) 



(Aa m -B)* =0, m/0, 



(16) 



where ^4 and -B are any operator not containing a+ or 
a_ m . In fact, for an appropriate operator C commuting 
with a m ,a_ m , a+,a± m , we have 

(Aa m B)y 

= (0| C*e e ™ a ™ a -™Aa m Be c ™ a ™ a -™C\0) 

pn k 

= E ^(0\C*(a m a_ m rAa m B(a+a±J k C\0) 



n,fc>0 



E 

n,k>0 



m 

n\k\ 



x (0\(a_ m nal m rC*ABC(a m r + \a+) k \0). 



Equation (|16[) follows because, on the one hand, 
{a m ) n+1 {a+) k \0) = if n + 1 > k and, on the other 
hand (0|(a_ m ) n (at m ) fe = if k > n. 

The total particle number is computed as 

( E a ™°™)* = ( a o a o )* + E ( a ™ a ™)* 



m#0 



° + E r 



m^O 



so the constraint (JSJ) is equivalent to 
TV = iV o + V 



(17) 



We shall see after (|5"o| that with our choice of parameters 
c m and -/Vq we have 



N- N 
N 



CX 3 ' 2 g 1 ' 2 + 0{\ 2 g)- 



(18) 



with a positive constant C that depends only on the 
unsealed potential Vq. In particular, the depletion rate 
of the condensate is of order \ 3 / 2 g 1 / 2 . 

We classify the interaction terms in the Hamiltonian 
(|3j) according to their number of zero momentum oper- 
ators, Oq or a . It will turn out that only even number 
of zero momentum operators give non-zero contribution. 
We will thus write 

V u a+a+a p _ u a q+u ^ = E + E 2 + E 4 , 

where Ek, k = 0, 2, 4, defined below, denote the contribu- 
tions of terms with exactly k zero momentum operators. 

Case 1. All four operators are with zero momentum, 
i.e. p = q = u = 0, and the contribution of this part is 



Eo = ^VoN 2 



(19) 



Case 2. By momentum conservation in the interaction 
term, it is impossible to have exactly three zero momen- 
tum operators. The terms containing exactly two zero 
momentum operators are 



1 1 u^O 



a o a u a -u + V u a+a± u a a 



(20) 



+ 2(K + V )a+a+a a % 



The contribution of this term to the potential energy of 
* is 



£ 2 = E 

p#0 



NoV p Cp -\- Cp 2(Vp + Vo)N \c p \ 2 
2|A| l-\c p \ 2 + 



2|A| 



\c P \ 2 



(21) 



Suppose that among the four momenta, p, q, p — u, 
q + u, exactly one is zero, say p — u (other cases 
are analogous). Then the remaining three operators 



5 



are ajala v+g . Since p,q,p + q are nonzero, either 



i p llg "p + g 

p,q,p + q, —p, —q, — (p + q) are all different, or p = q, 
in which case a p+q stands alone without any other oper- 
ator or a ±^ p+q y From (|16|) . the expectation of 
this term with respect to ^ vanishes. This proves that 
there is no contribution for the case of exactly one zero 
momentum operator. 

Case 3. No zero momentum operator is present in the 
interaction term in ([3]). Let r := p — u, s := q + u. 
Based upon only the following cases yield a non- 

zero contribution: 

• p = —q and r £ {p,~p}, but since r must be 
±s and we have momentum conservation, r = — s. 
The energy contribution is the main term in this 
case: 



p^O r^0,±p 



2|A|'E E Vp ~ r (l- - |< V | r i 



C75 Cf 



• p = = r = — s 

1 f> v-/ + + , _ ^0 \- kp| 2 (l + |cp| 2 ) 
^ pVo " " " " " 2|A| ^ (1-M 2 ) 2 



• p = —q — —r = s 



p#0 



-77 ^ ^2p( fl P a -p a ~p a p )* 



2|A| 



p#0 



^ V f> |c P | 2 (l + |cp| 2 ) 
2|A| ^ 2p (l-|cpP)2 



p = q, then r+s = 2p and r = ±s implies r = s = p 
and we have 



^ E 3><«a, «, >* - 2jxf E (T _ | Cp | 2)2 



2|cp| 4 



• p = 7', (7 = s and p — s, (7 = r but p 7^ ±q: 



p^O q^O, ±p 
1 

2IA 



1 1 p^Qq^0,±p 



p^O q#0 : ±p 



(l-| C p| 2 )(l-| C ,| 2 ) 



Collecting all these terms, the contribution of the case 3 
to the potential energy is 



E4 = E\ -\- f^4 



(22) 



We now combine the contribution to the potential en- 
ergy from case 1 and case 2 and we use the relation 
between ./V and Nq given by (fl7|) : 



Eq + Ei = 



(23) 



1 



~2|A| 



V N$ 



+E 

p#0 

1 

W\ 



N V P Recp 2(y p + y )iVo M 2 



|A| l-|cp| 2 2|A| 1-M 



^oiV 2 



v |Cm|2 

|A| 



— (y- 

2|A| V A'. 1 - |r 



2 s 2 



E 

p#0 

E 

p#0 



A% Recp (T/p + y )A^ |c p | ; 
|A| l-|c p | 2 + 



|A| l-|cp| 2 
Rec„ 



12 1 l-|cp| 2 



(V P + Vo) 
|A| 



Ei 



I Cm I 



1 - C 



Notice that there are two terms of the form 
~]T[^oX)m^o i-|cj |^ wri ich cancel each other. So we 



(24) 



can rewrite (|23j) as 

£-0 + £-2 = £() + £ 2 + "2 j 

where 

E = ^V N> , 



(25) 



^ 2 = E 

p#0 



f> M f> Rec, 



P i-|c P | 2 



l-|cp| 2 



I Cm I 



|A| V " J. - 



Re Cp 

1 - ici; 



(26) 



and ^2 is given in (|13p . The first term in (f!?6")) , when 
combined with the kinetic energy contribution 



p#0 



is the second term in (|15[) . The remaining main terms 
in (fTS"]) come from the rest of -E2, £0 and £4, i.e. 

Em = Ep 2 1 i C f 2 , 2 + £0 + £2 + £4 ■ (27) 
p#o 1 |c f | 



where was defined in (fl"4 



This completes the proof of Lemma [5J □ 
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Notice that the main terms in the potential energy 
come from the following channels: 

a+a+a a + ^ ^ (v u a+a+a u a_ u 

+ V u a+a± u a a Q + 2(V a + V )a+a+a a u 

+ 2U\T ^ H V p -rap a -p a r a -r ■ ( 2 8) 

The main energy contribution from these terms are all 
of order Ng. In the last term, the interaction between 
two large momenta (|p|,|r| ~ 1) pairs contribute with 
order Ng. The order Ng 3 / 2 term comes partly from 
substituting iVo, the expected value of a^a , with N 
(using that N — N ~ CNg 1 ^ 2 ) and partly from the 
interaction between a low momentum pair, |p| -C 1, and 
a large momentum pair, |r| ~ 1. 



THE ONE-PARTICLE SCATTERING PROBLEM 

Recall that 1 — w was the solution to the zero energy 
scattering equation ([6]) and w p denotes the Fourier trans- 
form of w. If V is smooth and it decays sufficiently fast 
at infinity, then w(x) is smooth and w(x) < C|a;| _1 for 
large \x\. Its Fourier transform on R 3 , J R3 e~ ipx w(x)dx, 
has a \p\~ 2 singularity at the origin. The lattice Fourier 
transform satisfies the regularized bound 

\w p \ < w = / w(x)dx < CL 2 , p e A* . 

J A 

This bound guarantees that for any function <p £ L 1 (R 3 ), 
the lattice Fourier transform of ipw can be computed as 

& w ) p = (<P * W)p = r^J ^ Vp-rW r 



In particular, from 



1 1 r^O 



Thus, modulo an error that is negligible in the thermody- 
namic limit, we can restrict the momentum summations 
involving w r to r ^ 0. 

From the definition of the scattering length (|7|), we 
have 



87Tfl = Vr 



' 



d 3 p ~ 

VnWr, — so = Vb - 



p p (2tt) 3 

From the scattering equation we get 
1 
2 



1 1 p#o 



1 

(29) 



2~ 1^ 1 - 

- p W p + -Vp - — ^ Vp- r W r 



Vp^O. 
(30) 



Define 

/(x):=^M4 g(x):=V(x)-f(x), 
then in Fourier space we have 

f p = (V* w) p = 7^7 ^-rWr + • 



(31) 



^7ra = 50 + O ( — 



(32) 



In the sequel we will not carry the negligible error term 
0(1/ L) in the formulas. 

Lemma 3 Let V(x) > 0, V ^ 0, be a radially symmetric 
smooth function with a sufficiently fast decay as \x\ — ► oo 
and let V = XV . Then, for a sufficiently small X, 

Vp,f p ,g p are real and have a fast decay as \p\ — > oo, 

(33) 



Moreover 



< f Q , go < Vo 



< A l = k = O(A) 
87ra g 



(34) 



(35) 



Furthermore, f p ,g p ,V p are uniformly Lipschitz continu- 
ous, i.e., 



sup (\g p - g p - r \ + \% - V^_r|J < CX\r\ ; 



sup|/p - f p - r \ < CX 2 



(36) 



with a constant C depending only on V . All statements 
hold uniformly in the thermodynamic limit, i.e. for all 
L sufficiently large. 

Proof. The reality of V p , f p ,g p follows from the radial 
symmetry. From the scattering equation (|30p 

-2p 2 w p +Vp-fp = 0, g p = 2p 2 w p Vp £ . (37) 

By iteration, we obtain the Born series for the scattering 
wave function (p =/= 0) 



Vp 1 x ^ Vp— r v r i x ^ 
V~V^ 2r 2 2p2 2^ 



Vp- r V r 1 



(38) 

It is easy to see from the expansion (|38p that (|3"3")) is satis- 
fied if Vp is sufficiently small and decaying. Furthermore, 

% = IAI-^^o^p = °( A2 )' while v o = cX with a 
positive constant c — J V, thus <?o = Vo — /o > cA/2 if 
A is sufficiently small and (|3"5|) follows. The rest of the 
statements of the Lemma also easily follows from 
□ 



THE MINIMIZATION 

From now on we assume that c p are real, it is most 
likely that complex choice does not lower the energy of 
our trial state. We introduce new variables 



r^0 



1 + C p 



7 



From \c p \ < 1 we have e p G (—00, |). We also have the 



p 

relations 

c 2 
°p 



1 



1 - 2e r , 



1 



l-2e„ 



(39) 



We shall choose e p via the following Lemma. This 
choice will become clear later on. 

Lemma 4 For any sufficiently small g < go{\) let 
—oo<e p <l/2be the minimizer of 



„2 C P 



1 - 2e„ 



eV P 



1 - 2e„ 



Then 



1-1 



V 2 + 2^/p 
and i/ie minimal value is given by 
1 r 



ofp^p 



1/2 



(40) 



(41) 



2 L 



(p 2 + 2^ p )(p 2 + gf p ) - (p 2 + g(V p + f p )) 



(42) 



Proof. Consider the minimization problem 

„2 



mm 

e<l/2 



1 - 2e 1 - 2e 



ce 



where the parameters satisfy a + 2c > and a + 2b > 
0. After differentiating in e, the minimizers satisfy the 
equation 



2(a + 2c)(e-e 2 ) + 6-c 



(l-2e) 2 

Solving the quadratic equation, we get 

b-c \ 1 / 2 







1 r 
2 



i± l 



a + 2c 



(43) 



and by the conditions on a, 6, c we have 1 + 2(6 — c)/(a + 
2c) > 0. In our case, since e < |, only the negative sign 
can be correct. With this choice, the minimum value is 



+b 



l-2e 



l-2e 
_ 1 

~ 2 



\/(a + 26) (a + 2c) - (a + 6 + c) 



In our application, the conditions a+2c > 0, a+2b > 
are equivalent to 



p 2 + 2gf p > 0, 



P 



2gV p > 



(44) 



and they are always satisfied if g is sufficiently small. In 
the regime \p\ > 4(f>t^j) 1 ' 2 , these inequalities follow from 
the bounds \V P \ < V and \f p \ < f < V (see pi])). For 
\p\ < 4(pVb) 1/2 wehavep 2 + 2g/ p = p 2 +2f?/ + O(e 3 / 2 ) > 
p 2 + gfo > by the Lipschitz continuity of f p (|36[) and 
the lower bound fo > 0. The other inequality in (|44f 
is proven analogously. Actually, these calculations also 



show that p 2 + 2gf p and p 2 + 2gV p have a positive lower 
bound uniformly in p, if g is sufficiently small: 

inf (p 2 + 2gf p ) > 0, inf (p 2 + 2g%) > . (45) 



□ 



We now rewrite the error terms f^2 + ^4 in terms of 



o 2 + n 4 

= E 

p#0 



e 2 e 2 



2|A| 

+ V 2p 



r^0,±p 

2e p J 
1 - 2e„ + 2e 2 



(1 - 2e p )(l - 2e r ) 



1 / e p (l-e p ) N 2 / 7 1-: 
All l-2e B J \ (1 



1 - 2e p + 4e 2 



(l-e P ) 2 



(l~e P ) 2 

% + Vo - 2 

|A| 



(y* e r \ e P 

\ 1 - 2e r / 1 - 2e r 



4- ^Wv^ e p 

2|A|V4-! l-2e. 



(46) 



p#0 



For the main term (fT5| . we replace Vo with 87ra + 
l A r X Ep^o ^p in ^0 (see (EH) and ([27])) by using © 
at the expense of a term of order 1/L that is negligible 
in the thermodynamic limit. Thus, neglecting this error 
term, we have 

g 2 - 
E M = AiraNg + — ^ V p w p 



p#0 



E 

p^O 



2|A| 4- P ~ r (1 



p ^"-p 
Cp(l Cp)e r (l e r ) 



r^0,±p 



Vp / ^ e r \ c p\ 

|A|V^l-2eJ~-2e 



(l-2ep)(l-2er) 



(47) 



By using (|34[) and (j3l)|) we have that for a sufficiently 
small but fixed S (depending on V), 

y < Vp < Vo, y < /p < Jo, y < ffp < So, 
for all |p| < <5. In particular we have 



<7o 



<7o 



(48) 



Using the lower bounds (|45[) and the approximation 
(|36|) of <?p, / p for small p, we obtain the following estimate 
on e p defined in (|4"Tj) : 



f Ie p | < C 



for Vp 

,M for ^VA 2 <N<5 (49) 



ego 



a I gl?pl 

-pi ^p 2 " 



for |p| > <5 , 
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where C is a constant depending on V and the notation 
A ~ B indicates that A and B have the same sign and 
\\A\ < \B\ < \A\. Here we have used the fact that < 

/o)<?o < Vo and that /o is order V 2 while g = Vo — fo- 
Similarly, we have 



and 



1 - 2e E 



Gp ( 1 6p ) 



< 



1 - 2 eiJ 



-TTW for W< 5 
for |p|>5 



for \p\ < 5 



(50) 



< 



p 2 +gfo 
el?pl 



(51) 



for H><5, 



where for positive quantities A < B indicates that A < 
CB with a constant C depending only on V. 

Lemma 5 Suppose e p is given by (|41[) . Then the energy 
E = (H)<j, of the state ^ satisfies 



E = AiraNg 

^-tV l-2e 



gV P 



V 



1 - 2e 



2|y\i E ^ _re P er + 9^ 

1 1 r^O 



P 

2 ^ 



(52) 



+ 0(Ng 2 \logg\) 

as g — > . 

Proof. We first prove that + f^4 are negligible. Note 
that, using \g p \ < g < Vo, and the bounds (|5D|) . 
we have 



-y- 



2e, 



|A| 



(£>5p) 2 



\p\>5 



< Cg*' 2 

Similarly, we find 

|ep(l-6p)| 



(53) 



1 \ - M 

IAI ^ 1 

1 1 p#o 



£.9o 



0^|p|<5P 2 + f?/0 



IAI V P 2 

1 1 H>« ^ 



(54) 



and 



lAI^I 1 
1 1 



e p (l - e p ) 2 < C 



2c r . 



Q9o 



— y (- 



< Cq- 



IAI — » 

1 1 \p\>s ' 

3/2 



(55) 



with constants depending on V. In terms of c p 's, we 
have, by (JM]), 



St 



?I*iiSC W .(M) 



p^O 1 Fl p^O 

The lower bounds in (|49| and (|50|l also imply that 
|^E 1 ^p-^ 3/ V /2 [l + 0(A,) 1 / 2] . (57) 



p#0 



In particular, we have proved (fT8|) after recalling (fT7| . 

Using (|53|) and ([55]) , the following terms are negligible 
from (146]) : 



TlH E (fc + ?p-r) 



2|A| 



p^O r^0,±p 



X 



e 2 e 2 



(1 - 2e p )(l - 2e r ) 



< CNg 2 



Vo 
2|A| 

E 



(Et^) 2 ^ 2 

p#0 



Up + Uo 
|A| 



(E i _ 2e r ) 1 - 2e„ ~ 



(58) 
(59) 

(60) 



\- (e P { 
where we used 



Vo f e p (l-e p ) ^ l-2e p + 4e 2 p f 



' (1 e p)^ 



<C^ a , (61) 



1 - 2e p + Ae 2 p 



= 1 + 3ct < 4 



(l-e p )2 



and similarly 



E 

p#0 



i> 2 p/ e p (l-e p ) N n-2e p + 2e 2 , 



|A| V 1 - 2e p J (1 - ep) 



< Cq 6 ' 1 . (62) 



Notice that the terms (|6Tj) and (|62|) are not extensive. 
All constants depend on V. We have thus proved that 
the f^2 + ^4 are bounded by Ng 2 . 

We now replace V p by V p ~ r in the last term of the 
main term Em j47j) . The difference can be estimated by 
using and (f5"l)) as 



Ew(EB-wa-)^ 



|A| 

p^O 1 1 r^O 



2e»r 



< 



C_ 

|A| 



E 



(^o) 2 |r| 



E 



{Q9r?\r\ 



^ s (r 2 + gfo)(r 2 + gV Q ) ££ 6 
\e p (l - e p )\ 

p/0 



xE^T 



26 E 



<CWe 2 |logf>| 



(63) 



We remark that this is the only term of size Ng 2 1 log g | 
and is the candidate for the third order term. 
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After this change, we can combine the two terms in 
the last two lines of (|4"7} as 

E l c% e p (l - e p )e r (l - e r ) 

2IAI ^ p ~ r fl - 2e„)(l - 2e r ) 

-ST — fSTv e r \ e P ( 1 - e p) 
^ |A| \^ p - r l-2eJ l-2e„ 



|A| 

p^O 1 1 r^O 



2|A| i 



(1 



e 

~2e 



_L V V — ^ 

2|A| p - r l-2e p l 

1 1 p.r^O P 



Cp ( 1 e p ) ^ ^ 



The last term comes from compensating for removing 
the restriction p ^ ±r. Similarly to (|61[) . it is not exten- 
sive and thus negligible. The second term is estimated 
exactly as (|59|) after estimating \V p - r \ < Vo; it is 0(Ng 2 ) 
and thus negligible. 
Since 

e p(l — e p) e p 

1 2c p 1 2c p 

this concludes the proof of Lemma [5l □ 

Proof of Theorem [lj Our goal is to minimize the 
energy given in (|52| in the parameters e p . This can 
be done either directly or via the following observation 
which converts the nonlocal term involving both e p and 
e r to a local term. For two functions ip, ?p defined on A*, 
we denote 

(V, V * ip) = ttj X) TaJ S Vp~ r ■ 

1 1 p^O 1 1 r^O 



Hence 



TXT 51 Taj 51 VrCper = (e, V * e) . 

We use the identity 

— (e, V * e) = — (e + £?u>, V * (e + gw)) — g(e, V * w) 



Q 



(w, V * w) 



and g — V — V * w (see (j3TjO . We can combine the last 
term with the last term in (1521) to obtain 



\Mq 2 



(jSs i v*€s)+^Y t ^ = \\MQ 2 (g,p~ 2 g) (64) 



p^O 



where we have used g p = 2p 2 w p . Thus we have 
E = inaNg 

p#0 



p 



1 - 2e„ 



2-"-2 



f?(^*w) p e p + 



Vr ^ P -r(e p + gw P )(e r + gw r ) 



2|A| 

oWliogg 



«2| 



(65) 



From the definition of e p in (|41|) . we have for p 2 > 



25P 



2(P 2 + 2g/ p ) 
Therefore, for p 2 > S _1 g, we have 

Q9p 



Hi) 



(66) 



e p + gw p 



2{p 2 + 2gf p ) 2p 



3 



= e 



dp + Iffp/p 
p 4 



(67) 



We now prove that the term in the last but one line 
of ((651) is negligible. By (J48J) and (gSJ we have 



2|A 



~n %-r(e P + QWp)(e 



p,r^0 
< 



C|A]( 



< CNg 2 



2 <S~ 1 g P 



Q9Q 
2 



gw r 



■ 2 >8- 1 g 



;5 P + \%fp\ 



The minimization of the first three terms in the sum- 
mation over p ^ in (|65|) is exactly given by Lemma 
[H With the choice (|4ip for e p , we have proved that the 
energy satisfies the following estimate: 



E = AirgNa + Q + 0(Ng 2 \\ogg\) 



(68) 



where 



Q 



P7 t0 L 



(? + ^)( ! 



2V, 



(69) 



We now take the limit L —> oo and change the sum- 
mation to integration. Changing variables, x = p 2 /g, 
d 3 p = 2ng 3 / 2 x 1 / 2 dx for x € (0, oo), then Q is given by 



irNg 3 / 2 f°° 
~(2^f 



F[x, y/xg)dx 



(70) 



where 
F(x,p) 



{x + 2j p ){x + 2V p )-(x+f p +Vp) + ^-)x 1/2 
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It should be noted that in our formulas, x and p are 
always related via the relation x = p 2 / g. 

With the notation a = 2f p /x, (3 = 2V p /x and recall- 
ing g p = V p — f p from (|57|) we can write 



F(x,p) = (V(l + a)(l+^) 

a + (3 (a-(3) 2 



r 3/2 



We divide the integration ([7(H) into x > eg 1 and x < 
cg^ 1 regimes for some small constant c. Since \V P \ + 
is bounded for all p, in the region x > cg^ 1 , we have 
|a| + |/3| -C 1 for c independent of g and g is small enough, 
a condition we assume from now on. We can thus expand 
a and /? in Taylor series and it turns out that the leading 
contribution is the third order term (a + (3)(a — /3) 2 /16. 
Hence we have F(x, yfgx) > for x > cg^ 1 and g small. 
Furthermore, we have the following estimate: 



/>oo />oo 

/ F{x, y/xg)dx < / x" 3/2 dx < Cy/g . 

Jcg^ 1 J eg -1 



(71) 



Similarly, the following inequalities, which will be useful 
later on, also hold: 

/>oo 

/ F{x, Q)dx < Cyftj, F(x,0)>0for x>cg- 1 . 

Jcg- 1 

(?2) 

For x < eg" 1 , we again use g p = V p — f p and rewrite 



F(x,p) 



G(x,p) 



-1/2 



G{x,p) where 



% [J (x + 2f p ) (x + 2V P ) + f p + V p -x] 



2[yf(x + 2f p )(x + 2V P ) + (x + f p + V p )] 
9 2 pW P + Vp)x-gl} 



2[^{x + 2f p )(x + 2V P ) + (x + f p + V p )] 
1 



(x + 2f p )(x + 2V P ) +x- fp-Vp 

The numerator in (|75|) may vanish only when 

x f p V p . 



(73) 



\x + 2f p ){x + 2Vp 
Solving for x, we have 



(U-v P ) 2 

4(/ P + %) 



(74) 



(75) 



In the regime x < eg" 1 , \p\ < s/c and from the continuity 
of f p , V p ([36]). the leading contribution of f p ,V p is given 
by /o, Vq > 0. Hence for c small enough (depending on 
A, but not on g), the solution (|75jl satisfies that 



x = {f * Vp T <fp + V p . 
^(f P + Vp) 



(76) 



Therefore, the numerator in ((73)) is positive in the neigh- 
borhood of the solution ([75]) and is thus also positive for 



all x < eg" . As a side remark, when combined with 
the previous argument for x > eg" 1 , this proved that 
G(x, yfgx) and F(x, yfgx) are positive everywhere. 

From (|73|) . G depends smoothly on x, f p and V p in 
the regime x < eg" 1 . Using the uniformly Lipschitz 
continuity of f p and V p ([36]) . we thus have 

\G(x, P )~G(x,0)\<C\p\(l + x)- 1 . 

Here we have used the second line of ([731 to obtain the 
decay in x for x large. Therefore, we have the error 
estimate 

pee' 1 peg- 1 i 

/ \F(x 1 ,/x^)-F(x,0)\dx<C v fg- / — — 
Jo Jo 1 + x 

< epilog g\ . 



Together with (J7TJ) and (f72|) . the same estimate holds if 
the integration domain is extended to the whole K + . 

From ©, (jnil> and (02) we see that h = %/g Q . With 
this notation we have 



( K \]{x + 2%){x + 2f + 2g ) 



( x + 2f +g ) + ^ 
2x 



+ 0(Ng 2 \\ogg\) 



x (y(y + 2h)(y + 2 + 2h) -(y + l + 2h) + 
+ 0(Ng 2 \logg\) . 



Recall go = 8wa from (|32|) and the definition of Q(h) 
from (fTTj) . We can thus write Q as 



Q = iiraNg 



^(h)(a 3 g)^ 2 

7T 



0(Ng 2 \\ogg\) , 



Together with (f6"5]) , this proves the main Theorem. □ 
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